Angular momentum projection of cranked Hartree-Fock states: Application to 

terminating bands in A~44 nuclei 
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We present the first systematic calculations based on the angular-momentum projection of cranked 
Slater determinants. We propose the /j,— >/ scheme, by which one projects the angular momentum 
/ from the ID cranked state constrained to the average spin projection of {ly) — I. Calculations 
performed for the rotational band in ^^Ti show that the AMP ly^I scheme offers a natural mech- 
anism for correcting the cranking moment of inertia at low-spins and shifting the terminating state 
up by ~2 MeV, in accordance with data. We also apply this scheme to high-spin states near the 
band termination in yl~44 nuclei, and compare results thereof with experimental data, shell-model 
calculations, and results of the approximate analytical symmetry-restoration method proposed pre- 
viously. 
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I. INTRODUCTION 

The energy density functional (EDF) method in nu- 
clear physics is nowadays the approach of choice for large- 
scale nuclear structure calculation. It has the same roots 
as the density functional theory (DFT) [1] in atomic 
and molecular physics, which is based on the Coulomb 
interaction between electrons. In nuclear physics, the 
EDF approach still lacks firm microscopic derivation and 
well-defined rules that would allow for systematic con- 
struction of an exact (or optimal) functional. In prac- 
tice, the nuclear EDF is constructed phenomenologi- 
cally, based on the knowledge accumulated within mod- 
ern self-consistent mean-field approaches built upon ef- 
fective density-dependent two-body interaction. 

These approaches, although successful in reproduc- 
ing gross nuclear properties and certain generic fea- 
tures of collective and single-particle nuclear motion are 
only coarse in confrontation with precise spectroscopic 
data. This situation calls for improvement, which can be 
achieved in two directions, namely, either explicitly, by 
using better and/or more complicated parameterizations 
of the nuclear EDF, or implicitly, by going beyond the 
mean- field or Hartee-Fock (HE) level. At present, intense 
studies are being conducted in both directions, aiming at 
exploring the limits thereof and answering the question of 
whether they can be considered equivalent, complemen- 
tary, or independent. In this work we will explore the sec- 
ond alternative, by employing the angular-momentum- 
projection (AMP) method of cranked HE (CHE) states. 

The spontaneous symmetry breaking (SSB) mecha- 
nism is inherently built into the HE approach. In many 
cases, it allows not only for incorporating a significant 
fraction of many-body correlations into a single HE state, 
but also serves as a source of deep physical intuition. 
Emergence of nuclear deformation (breaking of rotational 
invariance) leads naturally to the collective rotational 



motion, and is one of the most spectacular manifesta- 
tions of the SSB in nuclear physics. 

The CHE approximation treats collective (rotational) 
and intrinsic degrees of freedom on the same self- 
consistent footing. This fact is at the base of the success 
of this simple approximation, because in atomic nuclei 
both energy scales are strongly interwoven, and dramatic 
structural changes may take place along rotational bands. 
Terminating rotational bands [2] are the best examples of 
such possible changes, whereupon the collective rotation 
is followed by the total alignment of valence particles. 

Although the energies of rotational states are correctly 
reproduced by the CHE approach, and the changes of 
shape and pairing as well as the recoupling processes of 
individual nucleonic orbits are well captured, the price 
paid is high. Indeed, the resulting wave packets (de- 
formed CHE states) |$/ ) are well localized in the angular 
degrees of freedom and thus they are broadly spread over 
many angular momenta |<i>/ ) = '^j wi\I) , with only the 
average value of the projection of angular momentum on 
one the axes (the y axis in our case) being constrained. 
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Apart from strongly deformed states, this feature of 
the cranking approximation precludes applications of this 
formalism to compute transition rates, which constitute 
extremely valuable source of structural information. The 
demand for symmetry-restoration is therefore well mo- 
tivated. Starting from deformed CHE state | the 
goal can be achieved by projecting onto the eigenspaces 
of the angular momentum. Most of the calculations that 
have been performed so far invoke the angular momen- 
tum method applied to non-rotating states, see, e.g., 
Refs. [3-5] and the reviews in Refs. [6-9]. This limits 
their applicability to low spin states, where the influence 
of rotation on the intrinsic states can be neglected, and 
leads to an overestimate of the nuclear moment of inertia 
(Mol) as compared to the (realistic) cranking estimate. 
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After the ground-breaking studies in Refs. [10, 11], the 
AMP of CHF states has not been performed in mod- 
ern self-consistent calculations in nuclear structure. Re- 
cently, in Ref. [12], we presented results of such cal- 
culations for the test case of collective rotation of a 
well-deformed nucleus ^"''^'Gd. The AMP procedure we 
used has been implemented within the code HFODD 
[13, 14]. The calculational scheme proposed and tested 
in Ref. [12], dubbed hereafter ly^I scheme, assumes the 
AMP of spin component / of the self-consistent CHF 
state \^iy)-, which is constrained to the same mean value 
of the projection on the y axis, i.e., {ly) = I. It com- 
bines the simplicity of the self-consistent ID cranking 
approach, and its ability to reproduce the correct Mol, 
with the AMP after variation method. 

In this paper, within the same formalism, we present 
the first systematic calculations of rotational states along 
terminating bands in the A^AA mass region. The paper is 
organized as follows. Methods of calculation and results 
are presented in Sees. II and III, respectively. In particu- 
lar, the AMP of cranked states along the rotational band 
in ^^Ti is discussed in Sec. Ill A and the AMP of states 
near the band termination is presented in Sec. IIIB. Fi- 
nally, summary and discussion are given in Sec. IV. 



II. METHODS 



density matrices between rotated states: 
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The transition density matrix is also used for the density- 
dependent term. This is the only prescription available 
so far satisfying certain consistency criteria, formulated 
and thoroughly discussed in Refs. [5, 19]. 

Due to the overcompleteness of the \IMK) states, 
Eq. (4) is solved within the so-called collective subspace 
spanned by the natural states: 
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which are eigenstates of the norm matrix Nk'k having 
non-zero eigenvalues (n^ 7^ 0): 
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In practical applications, the cutoff parameter C is used 
and the collective subspace is constructed by using only 
those natural states that satisfy rim > C- By ordering 
indices m of the natural states in such a way that larger 
indices correspond to smaller norm eigenvalues, we can 
write the solutions of the HW equation (4) as: 



The angular-momentum-conserving wave function is 
obtained by employing the standard operator 
[15, 16] projecting onto angular momentum /, with pro- 
jections M and K along the laboratory and intrinsic z 
axes, respectively, 

\imk)=PIj,\^i^)='^-^ J D'^j^in) Rin)\^ijdn. 

(2) 

Here, fl represents the set of three Euler angles a(3^, 
(SI) are the Wigner functions [17], and R{0,) = 
yp-i-fiz jg ^]^g rotation operator. 
Since K is not a good quantum number, different 
K components must be mixed with the mixing coeffi- 
cients determined by the minimization of energy. The 
isT-mixing is realized in a standard way by assuming: 

17M)W =Y,9''^\IMK) ^ ^g'^^PLxm, (3) 
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and by solving the following Hill-Wheeler (HW) [18] 
equation: 



where Hk'k = mHPji'K\^) and TVj 
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denote the Hamiltonian and overlap kernels, respectively, 
and ^^'^ denotes a column of the g^^ coefficients. The 
overlap and Hamiltonian kernels have their standard 
functional form but depend upon transition (or mixed) 
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where the mixing coefficients of Eq. (3) read: 
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We can now define two types of the X-mixing. On 
the one hand, by the kinematic isT-mixing we understand 
the situation where only one collective state is used, i.e., 
"T-max = 1- Then, the solution of the HW equation 
amounts to calculating only one matrix element of the 
Hamiltonian kernel. 
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i-e-, f^^ = 1 and /, 

if- mixing, the mixing coefficients 17^' = ^/^VV^i ^'^^ 
entirely determined by the norm kernel and do not de- 
pend on the Hamiltonian kernel, i.e., they are entirely 
given by the cranking approximation and Coriolis cou- 
pling. On the other hand, by the dynamic X-mixing 
we understand the full solution of the HW equation for 
"T-max > 1, where the cutoff parameter Q is adjusted so 
as to obtain a plateau condition for the lowest eigenvalue 
El . Here, the generator-coordinate- method (GCM) mix- 
ing of different K components becomes effective, which 
potentially can modify the cranking mixing coefficients. 
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for m > 1. In the kinematic 
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We stress here that the kinematic _ft'-mixing docs corre- 
spond to a /f-mixcd solution too, and is not assuming 
any single given value of K. 

The deformed CHF states were provided by the code 
HFODD, which solves the Hartree-Fock equations that 
correspond to the Ritz variational principle, 



(11) 



with angular frequency u) adjusted so as to fulfill con- 
straint (1) and the value of ly being equal to /, accord- 
ing to our ly^I scheme. The y-signature and parity 
symmetries were conserved. The Hamiltonian and over- 
lap kernels were calculated using the Gauss-Chebyshev 
quadrature in the a and 7 directions and the Gauss- 
Legendre quadrature in the (3 direction [20]. In the nu- 
merical applications presented in this work we used the 
SLy4 [21] Skyrme force, but similar results were also ob- 
tained by using the SIII [22] force. The time-odd terms 
in the Skyrme functional were fixed by using values of 
the Landau parameters [23, 24]. The harmonic-oscillator 
basis was composed of 10 spherical shells. The integra- 
tion over the Euler angles was done by using a cube of 
50x50x50 integration points. 



III. RESULTS 




> 



0.4 
0.2 




1=2 4 6 8 1012 
A A»04 O 



Q Q 5 ^ O O O 0-0 ♦ 

I I I I I I X I I I 



9 m 



max 



In Ref. [12], we demonstrated that in a well-deformed 
nucleus ^^^Gd, the AMP /y— scheme reproduces well 
the cranking Mol along the rotational band, after taking 
into account the dynamic X-mixing. Here, we study the 
AMP methods applied to the terminating bands in the 
j4~44 region of nuclei. 



A. Ground-state rotational band in "'^Ti 

As a first example, consider the ground-state rotational 
band in "^^Ti. According to the CHF model, the shape of 
this nucleus undergoes a gradual change along the band. 
Starting from a well elongated {(32 ^ 0.23) shape at low 
spins, the nucleus goes through the alignment processes 
of the /2 protons and neutrons, and eventually reaches 
a nearly spherical {(32 ^ 0.05) shape at the terminating 
spin of / = 14, as shown in Fig. la. 

The calculated and experimental ^^Ti rotational bands 
are shown in Fig. lb. Results of the CHF calculations 
(left) are compared to those of AMP (middle) and to 
experimental data (right). In order to visualize the role 
of the JT-mixing, we have depicted results of the AMP 
calculations separately for the kinematic and dynamic 
if-mixing. For the dynamic X-mixing, the AMP excita- 
tion energies were calculated from the plateau condition. 
The stability of results with respect to the number of 
natural states mmax is shown in Fig. Ic. It can be seen 
that already at rrimax = 2, a perfect stability is obtained. 



FIG. 1: Results of the CHF and AMP calculations for the ro- 
tational band in '"'Ti obtained for the SLy4 interaction. The 
upper panel (a) shows the evolution of quadrupole and hex- 
adecapole shape parameters along the band as calculated us- 
ing CHF approximation. The middle panel (b) shows the ex- 
citation spectra calculated using the CHF approximation (left 
column) and AMP method with kinematic and dynamic K- 
mixing (middle columns), compared with the empirical data 
[25] (right column) . The lower panel (c) shows deviations AE 
of energy levels with spins I = 2, ... ,12 from the converged 
values corresponding the plateau condition, as functions of 
the number mmax of the norm eigenvalues used when solving 
the HW equation. 



i.e., here, the difference between the kinematic and dy- 
namic iC-mixing is related to adding the m — 2 state 
to the collective subspace. At higher values of nimax, 
one observes small departures from the converged val- 
ues, which are due to accidental mixing with spurious 
solutions (the largest such an effect is seen for / = 4 at 
Timax = 5). Nevertheless, physical converged solutions 
are clearly seen well beyond the point where the spuri- 
ous solutions become lower in energy, as is well known 
for other generator-coordinate- method calculations [26] . 

Figure lb clearly shows that within the ly^I scheme, 
the AMP effectively causes a decrease of the mean Mol 
within the band. This effect is expected to be generic 
for rotational bands of decreasing collectivity. Indeed, in 
such cases rotational correction is large at low spins and 
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FIG. 2: Similar as in Fig. lb, but shown in the absolute 
energy scale. Results of the AMP from the (ly) = ground 
state (first column) arc compared to those of the CHF (second 
column), AMP with kinematic and dynamic ii'- mixing (third 
and fourth columns), and CHFB calculations (fifth column). 



decreases at higher spins. This is illustrated in Fig. 2, 

where the calculated bands are shown in the absolute 
energy scale. In the case of ''^Ti, a net increase in the 
excitation energy of the terminating state, due to the 
AMP, amounts to about 2 MeV, but it is almost entirely 
related to the lowering of the 7 = state, while the 
terminating 7 = 14 state remains almost unaffected by 
the AMP. 

In the same Figure, we also show for comparison the 
AMP spectrum obtained by projecting from the {ly) = 
ground state. In this case, the overall Mol turns out to be 
much too small and the excitation energy much too high 



as compared to the CHF or AMP ly^I results. This 
fact once again shows the importance of the structural 
changes that occur in the system with increasing angular 
momentum. 

Note that the (ly) = state is axial, and thus for 
all spins contains only the K = component. There- 
fore, here, the collective space contains only one state 
and there is no TT-mixing at all. Moreover, axial sym- 
metry loads to a tremendous simplification of the AMP 
method, whereby only one-dimensional integration over 
the Euler /? angle is needed. Note also that the CHF 
states for (ly) ^ arc never axial, because the Cori- 
olis coupling always induces some non-zero nonaxiality. 
Therefore, the full three-dimensional integration over the 
Euler angles q;/37 is needed in our ly^I scheme. 

In the right column of Fig. 2, we also show results of the 
cranked Hartree-Fock-Bogolyubov calculations (CHFB) 
performed by using a zero-range volume-type (density- 
independent) interaction in the pairing channel. Its 



strengths of V„ 



-217.0 and V„ 



-237.5 MeV fm- 



for neutrons and protons, respectively, with the cutoff 
energy of Ccut = 50 MeV [27], was adjusted so as to re- 
produce, on average, the pairing gaps in this region of 
nuclei. For the (7^) = ground state of ^^Ti, the calcu- 
lated gaps read A„ = 1.394 and Ap = 1.632 MeV. 

One can sec that, in comparison with the CHF results, 
pairing correlations only affect states at low spins, 7 < 6, 



Pairing correlations vanish in the terminating state, and 
thus its energy is the same within the CHF and CHFB 
methods. At present, the code HFODD cannot perform 
the AMP of paired states, and we are yet unable to eval- 
uate combined effects of pairing and AMP from cranked 
states. 

It is also interesting to observe that, in the case of 
*^Ti, the dynamic TT-mixing is effective essentially only 
for 7 = 2 and 4, sec Fig. lb,c. This result is rather sur- 
prising, particularly in view of the fact that the only dis- 
continuity in spatial anisotropy of the CHF solutions can 
be seen around spins 7^6,8. This behavior is also qual- 
itatively different from that found in ^"''''Gd [12], where 
the dynamic TC-mixing was effective up to the highest 
(7 ~ 20) calculated spins. Apparently, the magnitude of 
TsT-mixing cannot be inferred solely from the shape but 
it also depends on individual (alignment) degrees of free- 
dom. 

In Fig. 3, we show probabilities Wj of finding the 
angular-momentum components in the intrinsic CHF 
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FIG. 3: Probabilities Wi (12) and Wk (13) of finding the 
given / and K components, respectively, in the intrinsic CHF 
states for different values of the projection on the y axis (ly). 
Solid lines are drawn only to guide the eye. Note that due to 
the conserved j/-signature symmetry, one has Wi=0 for odd 
values of I and W-k = Wk- 
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states, i.e., 

i'^h\P'KK\^Iy)- (12) 

K=-I 

The inset also shows probabilities Wk of finding the given 
K components, 

WK = {^I,\PkK\^I,)/Wi, (13) 

in one of the intrinsic states, {ly) = 6. At low spins, 
the Wi probability distributions follow the general pat- 
tern of deformed collective states. This pattern suddenly 
disappears at the termination point of (ly) = 14, where 
the CHF state becomes totally aligned, and therefore, it 
contains only the single / = 14 component. The Wk 
probabilities do not follow the standard pattern of col- 
lective rotation, where they correspond to the Coriolis 
mixing only, cf. Fig. 2 in Ref. [12]. Here, the aligning 
state causes the values of Wk to decrease with / and 
induces the crossing of those corresponding to even and 
odd K values in function of /. Note also, that the angu- 
lar momentum aligns along the y axis, while the standard 
projections K are calculated with respect to the z axis. 
Therefore, we do not expect to see a single K component 
of the aligned {ly) = 14 state. 

We conclude this section by analyzing excitation spec- 
tra calculated using the AMP of CHF states constrainc;d 
to different values of {ly), as shown in Fig. 4. It turns out 
that the AMP energies minimized with respect to {ly), 
or equivalently with respect to the angular frequency w, 
significantly differ from those obtained within our ly^I 
scheme. Indeed, although the minimum 7 = energy is 




FIG. 4: Excitation spectra calculated using the AMP of CHF 
states constrained to different values of {ly), with kinematic 
Jf-mixing. Solid lines connect states of given projected angu- 
lar momenta I and are drawn only to guide the eye. Arrows 
indicate the minimum energies in function of {ly)- For com- 
parison, the right column shows the CHF spectrum. 



obtained by projecting the {ly) = CHF state, this is 
not the case for / > states. For example, the mini- 
mum 1 = 2 and 14 energies are obtained by projecting 
the {ly) = and 10 CHF states, respectively. The result- 
ing excitation spectrum (in Fig. 4 indicated by arrows) 
is quite irregular and can hardly be associated with the 
physical result. The reason is the fact that the minimiza- 
tion over u> does not constitute the full variation-after- 
projection minimization of energy, which should be per- 
formed in the complete parameter space, including, e.g., 
the deformation parameters. Moreover, the fact that the 
I = \A state may gain energy by using collective correla- 
tions on top of the fully aligned (terminating) state points 
out a possible inadequacy of the EDF method when it is 
used within the variation-after-projection procedure, and 
not within the ly^I projection-after-variation scheme. 



B. Angular-momentum projection near the band 
termination 

In the vicinity of band termination, the number of con- 
tributing configurations drops down and the physics sim- 
plifies significantly. Reliable approximate analytical sym- 
metry restoration schemes can be easily derived for these 
cases; for details we refer the reader to the analysis pre- 
sented recently in Refs. [28, 29]. In the present paper we 
aim at further studying and testing these approximate 
methods against rigorous AMP results. 

Let us first consider the energy splittings between 
the favored- and unfavored- signature terminating states, 
[/7/2]^max and [/"/2]/max-i> rcspcctivcly, within the [fJf^i^] 
configurations, where n denotes the number of parti- 
cles in the /7/2 sub-shell. Within the nawe noncollec- 
tive cranking model, the unique aligned states 
can be considered as many-body reference states (HF 
vacua) with projections of the angular momentum be- 
ing conserved quantum numbers equal to the maximum 
allowed values, ly = /max- From these local HF vacua, 

the states can be generated by particle-hole 

(ph) excitations; in particular, by changing either the 
signature of a single neutron (r = v) or a single proton 
(r = tt). In spite of the fact that the underlying CHF 
solutions are almost spherical, they manifestly break the 
rotational invariance. Indeed, the two /max— 1 CHF solu- 
tions have conserved projections of the angular momen- 
tum, ly = /max— 1, but are in this case linear combinations 
of the total-angular momentum states with / = /max and 
/max — 1, i-e., up to a normalization factor: 

1^/ X — l) ^ ^|-^max; /max f ) ~t~ '^|-^max li /max 1)? 
\^f^ _l) a|/max;/max-l) - &|/max- l^/max-l)- 

The simplicity of the encountered situation allows for 
an approximate analytical estimate of mixing coefficients 
a and b [28, 29]. Equivalently, one can find these coeffi- 
cients by performing the exact AMP of the \^'f\ -i) '^^ 
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l^/lix-i'* states. The resulting probabilities P(/max — 1) 

of finding the /max— 1 components within the l^^J'^ _^) 
states are shown in Fig. 5a. The AMP results match 

perfectly the analytical results obtained in Refs. [28, 29], 
confirming reliability of the approximate method. 

Calculated energy differences E{Iynax) — E{I^ax ~ 1) 
are shown in Fig. 5b. Since the CHF solutions break 



the isobaric invariance, the AMPs of the |<I> 



and 



l^/max-i) states are not fully equivalent, and lead to 
slightly different energies. Results shown in Fig. 5 repre- 
sent arithmetic averages of both AMP energies. The only 
exception is ^^Sc, where we were able to perform numer- 
ical integration with a desired accuracy only when pro- 



jecting from the |$ 



CHF state, and the depicted 



point represents this single result. It is evident from the 
Figure that, except for ^^Sc and ''^Sc, the quality of the 
results is comparable to the state-of-the-art shell-model 
calculations of Refs. [29, 30]. 

The unfavored-signature /max— 1 CHF states discussed 
above were created by building the signature-inverting 
ph excitations on top of the /max reference state. Sim- 
ilar procedure applied to create /max — 2 CHF solutions 
leads to several nearly-spherical states, located, on aver- 
age, above the reference state. 

In this case, however, the SSB mechanism enters the 
game by pushing the collective (deformed) CHF solution 
down, below the reference state, and relatively close to 
the empirical energy. Hence, by going from the /max— 1 
to /max— 2 states, the physics changes quite dramatically, 
showing clearly two contrasting facets of the SSB mech- 
anism of rotational symmetry in nuclear physics. 

In the case of the /max^l states, the symmetry restora- 
tion results in a repulsion of two nearly- degenerate pro- 
ton and neutron CHF states, which are located above 
the reference /max state. On the one hand, the reorienta- 
tion mode, [/max; /max— 1), IS shifted down and becomes 
degenerate with the [/max; /max) solution, as required by 
the rotational invariance. On the other hand, the physi- 
cal mode, [/max — 1; /max — 1), is shifted up and becomes 
the unfavored-signature terminating state. 

In the case of the /max — 2 states, the SSB mechanism 
results in a repulsion of several nearly-degenerate pro- 
ton and neutron ph states. The collective CHF mode, 
l^/max-2), is shifted below the reference /max state, in ac- 
cordance with data. By the AMP mixing, the symmetry- 
restored collective mode projected from |$/„^x-2), i-C-, 
the 



1 — 2- 

-*max ^, 



/max — 2) state, gains some additional bind- 



ing energy. The situation described above is schemati- 
cally illustrated in the inset of Fig. 6. 

Calculated energy differences, /^(/max) — /^(/max — 2), 
are shown in Fig. 6. The CHF solutions, except for ^^Ca 
and ^^'^^Sc, correspond to collective states having /32 ~ 
0.10 - 0.12. The AMP shifts these states almost uni- 
formly down by about 300 - 400 keV, enlarging the split- 
ting by that amount, and improving an overall agreement 
between theory and experiment. It is, however, evident 
from the Figure that, here, the AMP does not improve 
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FIG. 5: Upper panel (a): probabilities of finding tlie /max — 1 
spin components in the CHF solutions. Squares 

and dots show results calculated using the exact AMP and 

approximate method of Rof. [28], respectively. Lower panel 
(b): energy differences between the favored- and unfavored- 
signature terminating states, [/7/2]-rmiix and [/7/2]-rmax-i! I's- 
spectively. Values calculated using the AMP method (di- 
amonds) are compared to those obtained within the shell- 
model [28] (circles) and empirical data [31-36] (dots). 



upon the incorrect isotopic/isotonic dependence of the 
CHF results. The magnitude of rotational correction is 
determined predominantly by the shape change, and does 
not vary from case to case. One can speculate that a de- 
tailed agreement with data would require an additional 
isospin-symmetry restoration. 

A similar trend was found for the /max — 2 states for 
configurations involving one-proton ph excitation from 
'^3/2 to /7/2 subshell, see Fig. 7. Here, all the lowest 
/max— 2 states are found to be deformed. The energy gain 
due to the AMP is again of the order of 300 - 400 keV 
and weakly depends on N and thus merely reflecting 
an increase of deformation between nearly-spherical /max 
states and terminating collective cranking solutions for 
/max — 2. This example confirms that the onset of col- 
lectivity causes a uniform mixing of various angular mo- 
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menta, regardless of individual features of specific nuclei. 
Indeed, in the studied nuclei, probabilities of finding the 
-^max— 2 spin component within the |$7^^^_2) = X)/ wi\I) 
CHF deformed wave packets are |w^/„ax-2p = 0.34±0.05, 
i.e., they extremely weakly depend on A'' and Z. More- 
over, values of \wi^_^^-2\'^ appear to be very similar to 
l'"'/maxP) which are equal to 0.35 ± 0.04, and again al- 
most do not change from one nucleus to another. 

Situation changes quite radically for the unfavored- 
signature [(^3/2/7/2^]^max-i states. Within the CHF ap- 
proximation, in iV Z and N Z + 1 nuclei there are 
three /max — 1 configurations that can be created from 
the /max reference state. Indeed, this can be done by 
a signature-inverting ph excitation involving either the 
neutron (j^) or proton (tt) /7/2 particle, or the proton 
^3/2 hole (tt), see Ref. [28]. The CHF solutions represent, 
therefore, mixtures of two physical |/max— 1; /max— i = 
1,2, states with the spurious reorientation |/max; /max— 1) 
mode. 

In such a case, the ly^I AMP scheme only removes 
the spurious mode, not affecting the mixing ratio of the 
two physical solutions |/max— 1; /max— i = 1,2. Hence, 
in contrast with the case of the [/7y2]^max-i states, qual- 
ity of the results strongly depends on the quality of the 
underlying CHF field. The AMP results corresponding 
to the 7f CHF solutions, which in '^'^-'^^Sc, 44-46rpj^ 
''^V are the lowest in energy, show that the admixtures 
of spurious components are of the order of 10%, see the 




44p„ 43Cp 45Cp 45TJ 46V 

*^Ca ''^Sc '''*Sc **Ti ^^Ti ^Ty 



FIG. 6: Energy difTcrcnccs between the [fy2]iniiix terminating 
states and the lowest [/7/2l/max-2 states. Dots and circles rep- 
resent the empirical data [31-37] and SM results of Ref. [28], 
respectively. Open diamonds label the lowest CHF solutions, 
which are collective except for the cases of *^Ca and *^'*^Sc. 
Full diamonds represent the AMP results. Relative energies 
between the /max reference state, non-collective ph /max — 1 
and /max— 2 excitations (thin lines), SSB effect in the /max— 2 
state (thick dashed line) , and final AMP configuration mixing 
(solid lines) are schematically illustrated in the inset. 



inset in Fig. 8. The obtained rotational corrections are, 
therefore, small of the order of 100 - 200 keV, and the 
disagreement with data remains quite large, as shown in 
Fig. 8. We show these results only as an example of possi- 
ble AMP calculations. However, for a complete analysis, 
one should, in principle, perform the GCM mixing of the 
AMP states corresponding to any possible CHF /max — 1 
configuration. A study in this direction is left for the 
future work. 



IV. SUMMARY AND DISCUSSION 

The present work reports on the first systematic cal- 
culations using the AMP of cranked Hartree-Fock states. 
The technique used, called the ly^I projection scheme, 
assumes the projection of the angular momentum compo- 
nent / from the one- dimensional (ID) cranked Hartree- 
Fock solution constrained to (ly) = I. The method ben- 
efits naturally from such nice physical features of the ID 
cranking model as the shape-spin self-consistency or the 
ability to get a realistic estimate of the nuclear Mol. It 
is shown that the ly^I AMP scheme leads to values of 
Mol that arc much more realistic than those obtained by 
using the AMP of non-rotating {ly) = states, as was 
the common practice up to now. 

In particular, application of the scheme to the rota- 
tional band in ^^Ti clearly improves the Mol at the bot- 
tom of the band. It also reveals a simple mechanism 
by which rotational corrections allow for improving upon 
an incorrect excitation energy of the terminating state, 
^-^/max ) obtained in the CHF calculations, which account 
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a 0.5 





42Ca 42Sc 44Sc 44Xi 46Xi il\ 
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FIG. 7: Energy differences between the [rf.j/2/7/2 l^i 
minating states and the lowest [cig^/j/ 



ter- 



-1 rri+l] 
3/2 J 7/2 



states. Dots 

and circles represent empirical data [31-37] and SM results 
of Ref. [28], respectively. Open and full diamonds label the 
collective CHF solutions and AMP results, respectively. 
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I I I I I I 

43SC 44Sc 45Sc 45Xi 46Xi ily 

FIG. 8: Energy differences between the [rf3/2/7/2^1^max termi- 
nating states and the lowest unfavored-signature terminating 
'^3/2/7/2^]^niax-i states. Dots represent empirical data [32- 
37] and open and full diamonds represent the CHF and AMP 
results, respectively. Inset shows probabilities of finding the 
^^max— 1 components within the lowest CHF solution | ^/^^^ - 1 ) . 

only for roughly half of the empirical value. 

Pairing correlations that arc active in the ground state, 
but do not affect fully aligned terminating state, can 
heal the situation only partially. Indeed, there is an 
upper limit for pairing correlation energy in the ground 
state, which can be sustained by a deformed A~44 sys- 
tem, equal to about 2 MeV. Further enhancement of pair- 
ing would induce the phase transition from deformed to 
spherical shape [38], and the rotational band could not 
have been built, contradicting the experimental data. 
Together, the rotational and pairing effects can bring 
^-^/max to about 9 MeV, i.e., some 10% below the 



experimental value. Let us mention here that for the 
/ = 14 state projected from the {ly) = groimd state 
one obtains Ai?/^^^~19 MeV, i.e., the result, which is 
well above the empirical excitation energy of the termi- 
nating state. 

For nearly-spherical unfavored-signature [/"/2]/max-i 
terminating states, our AMP calculations give results in 
excellent agreement with data, and validate approximate 
projection methods introduced in Rcfs. [28, 29]. We also 
show that the onset of collectivity in the [/7/2]/max-2 
states is quite correctly reproduced by the CHF calcu- 
lations, on top of which the AMP gives only a small cor- 
rection going in the right direction in comparison with 
data. However, details of the isotopic dependence are 
not reproduced here. 

Similar conclusions are obtained for the [<^3/2/7/V] ^o^' 
figurations that involve one-proton ph excitation across 
the Z = 2Q shell gap. In this case, in both /max — 1 and 
/max— 2 states near the band termination the collectivity 
sets in, while the energy differences with respect to ter- 
minating /max states are underestimated in the CHF and 
AMP calculations. 

The AMP of cranked Hartree-Fock states presented in 
this work was performed by applying the standard pro- 
jection techniques to the CHF solutions obtained within 
the EDF method. We have checked that all the results 
are stable with respect to numerical parameters such as, 
e.g., numbers of integration points used when integrating 
kernels over the Euler angles. Nonetheless, one should 
be aware of potential risks caused by difficult to control, 
uncompensated poles plaguing projection techniques of 
states obtained within the EDF method [12, 39]. Clearly, 
the future of projection methods crucially depends on a 
satisfactory solution of the problem of such singularities. 

This work was supported in part by the Polish Ministry 
of Science and by the Academy of Finland and University 
of Jyvaskyla within the FIDIPRO programme. 
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